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Abstract 

The letter proposes a procedure for generation of strongly chaotic beats that have been 
hardly obtainable hitherto. The beats are generated in a nonlinear optical system governing 
second-harmonic generation of light. The proposition is based on the concept of an optical coupler 
but can be easily adopted to other nonlinear systems and Chua's circuits. 
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I. INTRODUCTION 



In the last few years, an increase in the interest in the study of dynamical systems 
and design of nonlinear circuits generating chaotic beats has been observed [Grygiel & 
Szlachetka, 2002; Cafagna & Grassi, 2004-2006; Sliwa et. al, 2007]. Beats in a nonlinear 
oscillator appear if the oscillator becomes quasi-periodic for any reason, for example, due to 
the frequency detuning [Minorski, 1962]. Usually, in nonlinear systems, we observe repeated 
sequences of complicated revivals and collapses [Eberly et. al, 1980; Abdalla et. al, 2005]. 
Chaotic revivals and collapses, in contradistinction to the quasi-periodic ones, do not have 
any repeated structure. The degree of chaoticity of beats generated in a dynamical system 
is usually confirmed with the Lyapunov exponents. This approach seems to be necessary as 
intricate quasi-periodic beats frequently resemble chaotic beats and vice versa. In particular, 
the numerical procedure proposed by Wolf [Wolf et. al, 1985] is a useful and efficient method 
to get such exponents (ordered from maximal to minimal value), known also as the spectrum 
of Lyapunov exponents. The first exponent is traditionally named the maximal Lyapunov 
exponent (MLE), and its positive value is indicated of a chaotic motion. Two positive 
Lyapunov exponents denote hyperchaotic behaviour. One way to get chaotic beats from 
their quasi-periodic counterpart is to decrease damping in the quasi-periodic system. This 
action, usually, leads to generation of weakly chaotic beats [Sliwa et. al, 2007]. Therefore, 
the problem is to find a way of generation of strongly chaotic beats. In this letter, we 
propose a solution to this problem, by way of example, in a simple nonlinear optical process 
- second-harmonic generation of light, using the concept of a typical optical coupler [Pefina, 
Jr. & Pefina , 2000] that is made of two identical systems interacting linearly with each 
other. Let us first consider the generation of chaotic beats in a single system. 



II. BEATS IN SINGLE SHG-SYSTEM 



We study the following dynamical system [Drummond et. al, 1980; Mandel & Erneux , 
1982; Gao, 2004]: 

^ = -iuja - 7a a + ea*b + Fe~ iQat , (1) 
^ = -i2ub - 756 - 0.5ea 2 , (2) 
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where the complex variables a and b represent the amplitudes of the fundamental and 
second-harmonics modes, respectively. The parameter e governs the nonlinear interactions 
between both modes. The quantities u and 2uj are the frequencies of the fundamental and 
second-harmonic modes, respectively. The terms 7 a a and 7 6 describe the mechanism of loss. 
Moreover, the system is pumped by an external field Fe~ lQai , where F is an electric field 
amplitude at the frequency fl a - The parameters, u>, j a , 7&, e, F, and fl a are taken to be real. 
If 7 a >> 7b, the loss mechanism in the second harmonic mode is usually neglected, that 
is we assume 7 = 0. This simplification is frequently called a good frequency conversion 
limit. The system ((T])-(j2J) has three periodic solutions. Namely, two coexisting solutions, 
when 7h = and Q± — uj ± 0.5F/7 a [Sliwa et. al, 2007]: 

a±(t) = -e-^\ (3) 

7a 

b±(t) = T'-e-^, (4) 

2 7a 

and the resonance solution, when Q = u [Mandel & Erneux , 1982]: 

a(t) = (A + B)e~ iuJt , (5) 

b(t) = -^{A + B) 2 e- i2u}t , (6) 
27fe 
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We may expect that revivals and collapses (beats) appear in the neighbourhood of the 
periodic solutions if the conditions of periodicity are not held. For solutions (J3j)-(jll) the 
region of beat frequencies is situated between the frequencies f2_ and fl + (see Fig. 2 in [Sliwa 
et. al, 2007]). Beats in this region are either quasi-periodic or chaotic depending on the 
value of the damping constant 7 a . Beats can also be generated in the neighbourhood of 
solutions ©-([6]) if Q a ~ uj. However, the generation of beats when the second harmonic 
mode is damped (7^ 7^ 0), is much less efficient than in the good frequency conversion limit, 
where 7^ = 0. Let us shortly consider this problem by way of a numerical example. The 
results are presented in Fig.l. The Lyapunov map shows the regions of quasi-periodic beats, 
chaotic beats and purely periodic behaviour; individual colours correspond to the values of 
the maximal Lyapunov exponents. The case of beats generation corresponds the inside of 



the parabola , while the periodic states correspond to the region outside it. As seen, the 
regions corresponding to chaotic beats (positive value of MLE) are in the lower part of the 
parabola ( blue and yellow areas). MLE's are of the rank 10~ 3 . Therefore, beats generated in 
the neighbourhood of the resonance (Q a = uj) are weakly chaotic. An attempt at enhancing 
the chaoticity of the beats by increasing the amplitude F has proved ineffective. Assuming 
e.g. F = 20 instead of F = 5 the chaoticity increases but the structure of the beats is 
destroyed. Therefore, the problem arises how to force the system (a, b) to generate strongly 
chaotic beats in the presence of large damping that allows only generation of periodic states 
or quasi-periodic beats. In what follows, we propose a solution to this problem. 



III. BEATS IN SHG-COUPLER SYSTEM 



Let us consider the system (a, b) and its copy (A, B) and join both systems linearly in 
the following way: 

= -iua - 7 a a + ea*b + Fe~ in<lt + s x A , (7) 

^ = -i2ujb - j b b - 0.5ea 2 + s x B , (8) 
dt 

dA 

— = -iuA--f A A + eA*B + Fe- iQAt -s 2 a, (9) 
d R 

— = -i2uB --f B B -0.5eA 2 - s 2 b. (10) 

(JjL 

In nonlinear optics, systems of this type are usually called two-core couplers. Therefore, 
the subsystem (a, b) interacts with the subsystems (A, B) via the si-terms and vice versa 
(A, B) interacts with (a, b) via the S2-terms. If, for example, S2-terms are turned off then 
the interaction is unidirectional and, consequently, (a, b) plays the role only of a receiver 
whereas (A, B) is simply a transmitter. This type of interaction can be useful to generate 
chaotic beats within (a, b) even if the free system (a, b) is not able to generate (due to the 
large damping) chaotic beats itself. Let us consider two simple examples. 



A. Generation of chaotic beats in a periodic system forced by external quasi- 
periodic beats 

Let us first consider the dynamics of the system described by equations (!7])- (jT0l) if the 
interaction between the subsystems (a,b) and (A,B) is turned off (si = s 2 — 0). Then, the 
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system (a, b) for uj = 10, j a = 0.5, 75 = 0.005, e = 0.1, F = 5, fi a = 10.5 and for the initial 
conditions a(0) = 10, b(0) = —5i, generates a periodic state (or tends to this periodic state, 
as a steady state, if it starts from arbitrary initial conditions). The periodic behaviour in 
Re 6(t)-component is shown in Fig. 2a. 

The system (A, B) for u = 10, 7^ = 0.5, 7 B = 0.05, e = 0.1, F = 5 and Q A = 9.9 
and for the initial conditions a(0) = 10 and 6(0) = —5i, generates quasi-periodic beats (the 
spectrum of Lyapunov exponents has the form {—0.0396, —0.2203, —0.2214, —0.2346}). The 
beats are presented in Fig. 2b. 

To get chaotic beats, we now turn on the coupling in ©-(JED unidirectionally , in such 
a way that the subsystem (a, b) plays the role of a receiver (si 7^ 0) whereas (A, B) is a 
transmitter (S2 = 0). Then, the quasi-periodic beats are transmitted from the subsystem 
(A, B) into the periodic subsystem (a, b) but not vice versa. The coupling can be turned on at 
an arbitrarily requested time t. For t = 40, si = 1.8 and S2 = the initially periodic system 
(a, b) begins to generate chaotic beats (Fig. 3c). The choticity of the beats is confirmed by 
the spectrum of Lyapunov exponents {0.0436, —0.1311, —0.4905, —0.8791}. Therefore, the 
MLE (Ai = 0.0257) is two orders of magnitude greater than its counterpart obtained in 
Section II. 

The situation presented in Fig.2c changes rapidly if at the time t = 40, the feedback 
(s 2 7^ 0) is turned on between the subsystems (a, b) and (A,B). Numerically, we have put 
si = S2 = 1.8. The new geometric structure of beats in Re 6(t)-component is now presented 
in Fig. 2d. As seen, after t > 95 the transient behaviour vanishes and finally we observe 
beats. Also the structure of the beats in the subsystem (A, B) changes due to the mutual 
interactions between the subsystems (a, b) and (A, B) . After t > 90 we observe in ReB(t)- 
component repeated sequences of two different beats (Fig.2e). It is obvious that beats in 
Figs.2d-e are typically quasi-periodic, in contrast to those in Fig. 2c. The lack of chaoticity is 
confirmed by the spectrum of Lyapunov exponents {—0.1555, —0.1556, —0.1559, —0.1560}. 
Therefore, the feedback interaction between (a, b) and (A, B) damages the chaos within 
beats in the subsystem (a,b). 

The degree of stability of beats for different values of the coupling parameter si if S2 = 
is shown in Fig. 3 (blue line). Beats appear for s± > 0.1. For 0.1 < s% < 1.43 the beats 
are quasi-periodic, and in the range 1.43 < si < 5 they behave chaotically. In this region 
the maximal Lyapunov exponent Ai increases exponentially with increasing value of the 
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coupling parameter si, and generated beats become more and more chaotic. The numerical 
investigation shows that the system generates chaotic beats since si < 10. For si > 10 the 
beats disappear at all but the system still remains chaotic. 

It is worth noting that all the maximal Lyapunov exponents in Fig. 3 become negative if 
we turn on additionally the feedback interaction (s 2 = «i) between the receiver and the 
transmitter. Both subsystems still generate beats but they are not chaotic. 

B. Generation of chaotic beats in a quasi-periodic system forced by other quasi- 
periodic beats 

The periodic behavior presented in Fig. 2a does not appear if, instead of the external 
frequency Q a = 10.5, we put Q a = 9.84. Physically, it means that the detuning between 
the frequencies uj = 10 and Q a is now much smaller. And consequently, in the subsystem 
(a, b) instead of the periodic behavior, beats appear (Fig.4a). Let us now suppose that the 
subsystem (A, B) is in the same state as in Subsection A (Fig.2b) and try to generate chaotic 
beats in the unidirectional regime (A, B) =>- (a, b). For t = 40, Si = 1.8 and s 2 = the 
receiver (a, b) generates beats having complex geometrical structure but they are not chaotic 
which is confirmed by the Lyapunov exponents {—0.0221, —0.1550, —0.4972, —0.7829}. The 
coupling is simply too week to generate chaos within beats. However, a stronger cou- 
pling allow us to generate chaos. For example, for s-y = 2.75 and s 2 = we observe 
distinctly chaotic beats in the subsystem (a, b) which is confirmed by the Lyapunov coef- 
ficients {0.0257, -0.0867, -0.5607, -0.8354}. As follows, the MLE (A x = 0.0257) is of the 
same order as that specified in Subsection A. The geometric structure of beats in Re b(t)- 
component is presented in Fig. 4 .The degree of chaoticity of beats for different values of 
the coupling parameter s\ (if s 2 = 0) in the range < S\ < 5 is illustrated in Fig. 3 (green 
line). Quasi-periodic beats occur for < s\ < 2.62 and 2.81 < s\ < 2.92. Chaotic beats are 
generated in the range 2.63 < s\ < 2.81 and 2.97 < s\ < 5. Similarly as in Subsection A, 
on turning on the feedback interaction (s 2 — S\ — 2.75) both subsystems generate beats but 
they lose their chaotic nature. 
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IV. FINAL REMARKS 



We have proposed a procedure to obtain beats of a higher degree of chaoticity than that 
hitherto obtained. Let us emphasize that the method presented (unidirectional coupling 
between a receiver and a transmiter) can be easily realized in arbitrarily dynamical systems. 
The degree of chaos in the beats generated in the proposed optical transmiter-receiver system 
is of the same order as in the Rossler system [Rossler, 1979]. The coupling transmiter-receiver 
could be easily applied also in Chua's circuits, where weakly chaotic beats have also been 
obtained [Cafagna k Grassi, 2004-2006]. 
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V. FIGURE CAPTIONS 



Figure 1. The values of the maximal Lyapunov exponents marked by separate colors for 
the system ©-© with u = 10, e = 0.1, 7 a = 0.5, F = 5, < -f b < 0.2 and 9.5 < Q a < 10.5. 

Figure 2. Time evolution of Reb(t) and ReB(t) for the system (171)- (11 01) . where uj = 10, 
7a = 0.5, 7 6 = 0.005, tt a = 10.5, 7a = 0.5, 7s = 0.05, tt A = 9.9, e = 0.1, F = 5, and: 
(a)-(b) Si = s 2 = 0; 

(c) si = 0, s 2 = if t < 40 and s 1 = 1.8, s 2 = if t > 40; 

(d) -(e) Si = s 2 = if t < 40 and s x = s 2 = 1.8 if t > 40. 

The initial conditions are: a(0) = 10, 6(0) = — 5i, A(0) = 10 and B(0) = —5i. 

Figure 3. The values of maximal Lyapunov exponents Ai, for the system fl71)- f|TUl) as a 
function of the coupling constant si, where s 2 = 0, uj = 10, 7a = 0.5, 7 & = 0.005, 7 a = 0.5, 
7b = 0.05, e = 0.1, F = 5, tl A = 9.9 and Q a = 10.5 (blue line), tt a = 9.84 (green line). 

Figure 4. Time evolution of Reb(t) and ReB(t) for the system (JTj)- (jl0p . where uj = 10, 
7a = 0.5, 7 6 = 0.005, Q a = 9.84, 7a = 0.5, 7 B = 0.05, Q A = 9.9, e = 0.1, F = 5, and: 

(a) si = s 2 = 0; 

(b) s 2 = 0, si = if t < 40 and s 1 = 2.75, s 2 = if t > 40. 

The initial conditions are: a(0) = 10, 6(0) = -5z, A(0) = 10 and 5(0) = -hi. 
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FIG. 4: 
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